00 

o 

O 
(N 



X 



Minimization variational principles for acoustics, 

elastodynamics, and electromagnet ism in lossy 

inhomogeneous bodies at fixed frequency 



Graeme W. Milton 

OO ' Department of Mathematics, University of Utah, Salt Lake City UT 84112, USA 

Q_^. Pierre Seppecher and Guy Bouchitte 

^ ■ Institut de Mathematiques de Toulon, 

?M ' Universite de Toulon et du Var, BP 132-83957 La Garde Cedex, France 



K^ . Abstract 

2^ ■ The classical energy minimization principles of Dirichlet and Thompson are 

^r^ , extended as minimization principles to acoustics, elastodynamics and electromag- 

netism in lossy inhomogeneous bodies at fixed frequency. This is done by building 
t~^ ■ upon ideas of Cherkaev and Gibiansky, who derived minimization variational prin- 

^^ , ciples for quasistatics. In the absence of free current the primary electromagnetic 

f— >) I minimization variational principles have a minimum which is the time-averaged 

electrical power dissipated in the body. The variational principles provide con- 
straints on the boundary values of the fields when the moduli are known. Con- 
versely, when the boundary values of the fields have been measured, then they 
C^ ' provide information about the values of the moduli within the body. This should 

have application to electromagnetic tomography. We also derive saddle point vari- 
ational principles which correspond to variational principles of Gurtin, Willis, and 
Borcea. 
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1 Introduction 

As our goal is to extend the classical minimization principles to acoustics, elastodynam- 
ics, and electromagnetism, let us begin by recalling the well-known classical results of 
electrostatics. Consider the equation of electrostatics 

V-eVr = 0, (1.1) 



for the real valued potential V^(x) in a body F containing a locally isotropic material 
with scalar real positive dielectric constant ^(x). This equation can be rewritten in the 
equivalent form 

D = eE, E = -VV, V • D = 0, (1.2) 

where the first relation is known as the constitutive relation linking the electric displace- 
ment field D(x) and the electrical field E(x) where these two fields satisfy the differential 
constraints embodied in the second two relations. Two variational principles are well 
known. We have the Dirichlet variational principle that 

W{V) = inf W^(i:) where W{V) = /e|VZr/2, (1-3) 

y- Jt 

in which the infimum is taken over all trial potentials V(x) satisfying the boundary con- 
dition that li(x) = V^(x) on the boundary dV. We also have dual Thompson variational 
principle that 

W{B) = inf WCD) where W(B) = [ |D|V(2e), (1.4) 

D Jr 

V-D = 

where the infimum is taken over all divergence free trial displacement fields D satisfying 

the boundary condition that D ■ n = D ■ n on dT. The quantity e| VV^p/2 = |Dp/(2£) 

is the electrostatic energy density inside the body, and so these variational principles 

are also known as energy minimization principles. In the mathematically analogous DC 

electrical conductivity problem, where e(x) and D(x) are replaced by the conductivity 

o'(x) and the current J(x), this quantity is half the electrical power dissipation and in 

that context the variational principles could more accurately be called power dissipation 

minimization principles. 

An appealing feature of these variational principles is that the minimum values W^(y) 

and W^(D) can be expressed in terms of the values of V{'x.) and D ■ n at the boundary 

dT: 

W{V) = W(D) = -- f (D- n)\/(x). (1.5) 

2 Jar 

So these variational principles provide a constraint on the values of Vo(x) = V{x.) and 

go(x) = — D(x) ■ n at the boundary dT. The relationship between these boundary 

functions go and Vq is linear and the relationship go = NVq defines the Dirichlet to 

Neumann map N, which is a real self-adjoint positive semi-definite operator. Therefore 

N is determined by knowledge of the quadratic form 

w{Vo) = {Vo,NVo)/2, (1.6) 

where for any two scalar or vector, real or complex, valued functions p(x) and q(x) 
defined on the boundary dT, we denote 

(p,q)= / p(x)-q(x). (1.7) 

Jar 



Given trial potentials Vl(x, Vq), parameterized by the function Vq, with V_{yi, Vq) = Vo(x) 
on dr, the Dirichlet variational principle implies the inequality 

w{Vo) < W{V), (1.8) 

and provides an "upper bound" on the map N. By "upper bound" we mean an upper 
bound on the associated non- negative quadratic form w{Vo). 

To obtain a "lower bound" one uses the Thompson variational principle. Since 
w{Vo) = when Vq is constant the null-space of N consists of all constant fields and 
so N is determined by knowledge of the quadratic form 

w{qo) = {qo,^-'qo)/2, (1.9) 

for all fields go(x) with zero average value on dT. The functionals w{Vo) and w{qo) are 
Legendre transforms of each other, 

w{Vo)=mi[{Vo,q)-w{q)], (1.10) 

q 

where the infimum is over all fields g(x) with zero average value on dT, and so an upper 
bound on the functional w{q) provides a lower bound on w{Vo). Given trial displacement 
fields D(x, go) parameterized by the function q^, with — D(x) • n = go(x) on dT, the 
Thompson variational principle implies the inequality 

w{qo)<W{B), (1.11) 

and provides a "lower bound" on the map N. 

Thus the Dirichlet and Thompson variational principles allow one to bound the Dirich- 
let to Neuman map N. As shown by ?) they can also be used in an inverse fashion: if N 
and hence the functionals w{Vo) and w(go) are known one can obtain information about 
the permittivity distribution £:(x) from the inequalities (11.81) and (II. lip by appropriately 
chosing the trial potentials li(x, Vq) and trial displacement fields D(x, go)- It has been 
proved that knowledge of N uniquely determines e{x.) when e{x.) is piecewise analytic 
(?), or more generally when ^(x) is smooth (?). 

In this paper we extend these variational principles to acoustics, elastodynamics, and 
electromagnetism, in lossy inhomogeneous bodies at fixed frequency. By lossy we mean 
that the time-averaged dissipation of mechanical or electrical energy into heat is pos- 
itive everywhere, which means that the imaginary part of certain tensors, such as the 
permittivity tensor, are positive definite. Like the Dirichlet and Thompson variational 
principles, the variational principles we obtain are minimization (not saddle point) vari- 
ational principles. Other variational principles which extend the Dirichlet and Thomp- 
son variational principles to dynamics, not assuming time harmonicity, were derived for 
elastodynamics and electromagnetism by ?) and Willis (?, ?). These are minimizing 
variational principles in the Laplace domain but not in the frequency domain. We will 
recover the Gurtin- Willis variational principles in the time harmonic setting where they 



correspond to saddle point variational principles. We remark that ?) and Willis (?, ?) 
derived variational principles which are dynamic analogs of the Hashin-Shtrikman (?, ?, 
?, ?) variational principles, involving trial polarizations fields. We do not consider such 
variational principles here. 

Of course given a (inhomogeneous) linear differential operator A, mapping ?Ti-component 
complex vector fields to m-component complex vector fields, there is a trivial minimiza- 
tion variational principle associated with the equation ^u = h for the rra-component 
potential u, in which h is a rri-component source term. When the equation .4.U = h with 
u = Uq on dr has a unique solution for u, the infimum in 



= inf / (^u - h) ■ (^u - h), (1.12) 

ii ir 

where the overline denotes complex conjugation and the infimum is over all fields u 
with u = Uq on c?r, is clearly attained when u = u. What separates such variational 
principles from those of Dirichlet and Thompson is that, as explained above, the latter 
variational principles give us useful information on the boundary values of the fields. We 
desire variational principles with this property. Also, from a mathematical viewpoint, 
the existance of a solution cannot be deduced from (11.121) . and in case there is a solution 
no interesting inequality results from (I1.12p . 

In an important development ?) extended the Dirichlet and Thompson variational 
principles to the quasistatic electromagnetic equations (and also to the quasistatic elas- 
todynamic equations) where the equations (II. ip and (11.20 still hold but ^(x), V^(x), E(x) 
and D(x) are complex, and the imaginary part of £:(x) is positive. Since we build upon 
their work let us briefiy review their ideas in this context. 

By rewriting the constitutive law D = eE in terms of its real and imaginary parts 

s:)K"-:0(-)' '-' 

where the prime denotes the real part, while the double prime denotes the imaginary 
part, it is evident that since ^"(x) > one has the Cherkaev-Gibiansky saddle point 
variational principle 

Q{V', V") = inf supQ(]/',K") = supinf Q(]/',]/"), (1.14) 

V' y't y" V' 

where 

«E',r)^X^(^i:).(::::„)(^^,). (..5) 

and the infimum and supremum are over trial potentials with V2'(x) = V^'(^) ^^^ Y^'i^) — 
V'(x) on the boundary dT. To obtain a minimization variational principle from this 
saddle point variational principle ?) make a partial Legendre transform of the saddle 



shaped quadratic form to convert it into a convex quadratic form. This is equivalent to 
rewriting the constitutive law in the form 

^E")=K-D')' <"" 

where the matrix 

s" + {sr/{e") e'/e"\ 

e'/e" l/e" ) ^ ' 

is positive definite when £:"(x) is positive. From this it is evident that one has the 
Cherkaev-Gibiansky minimization variational principle: 

Y{V',B') = mf inf F(Z',D'), (1.18) 

V ■ D' = 

where 

nr,»^/(f')..(f'), (1.1.) 

and the infimums are over all trial real potentials with V^{^) = V^'(x) on dV, and all trial 
real divergence free displacement fields with D ■ n = D ■ n on SF. 

As an aside we remark that, as shown in section 18 of ?), one can extend these ideas 
of Gibiansky and Cherkaev to obtain variational principles for problems where the consti- 
tutive law takes the form D = eE, where the matrix valued field £(x) is not symmetric. 
For example, for conduction in a fixed magnetic field, the conductivity tensor cr(x) is 
real but not symmetric. Expressing cr(x) = crs(x) + cra(x), where crs(x) = crj'(x) and 
<Ta(x) = — cr^(x) are the symmetric and antisymmetric parts of the conductivity tensor, 
one considers current fields J(x) and electric fields E(x) which solve the conductivity 
equations 

J = (cr, + 0-JE, E = -\/V, V-J = 0, (1.20) 

in conjunction with current fields J(x) and electric fields E(x) which solve the adjoint 
problem 

J = (cr, - cr„)E, E = -VV, V ■ J = 0. (1.21) 

By adding and subtracting the two constitutive laws one gets a new constitutive law 



Js \ / CTs CTa \ / E^ 



-3a / \ — cr„ — cr, / \ E^ 



:i.22) 



involving fields 

J, = J + J, Ja = J-J, E, = E + E, Ea = E-E, (1.23) 

that satisfy the differential constraints 

V-J, = V-J, = 0, E, = -VK, E, = -VK, (1.24) 



where Vs = V + V and Va = V — V. Conversely, given any fields which solve f ll.22p 
and (ll.24p . one can recover the associated fields which solve f ll.20p and fll.2ip . The 
two problems are equivalent. Since the matrix entering the constitutive law fll.22p is 
real and symmetric there is clearly an associated saddle point variational principle when 
^^(x) is positive definite for all x. ?) show how this variational principle provides useful 
asymptotic estimates for the mathematically analogous problem of convection enhanced 
diffusion. Also by rewriting the constitutive law (11.221) in the equivalent form 






;i.25) 



which involves a matrix which is symmetric and positive definite when crg(x) is positive 
definite for all x, one obtains an associated minimization variational principle (?). In 
this paper we do not treat such cases where the constitutive laws involve non-symmetric 
matrices: this will be covered elsewhere. 



2 General Theory 

Let d be the dimension of the space and let J^(x) and ^(x) be two complex valued fields 
of the form 

where F(x) and G(x) are m x d dimensional matrix fields and f(x) and g(x) are m 
component vector fields, and which satisfy the differential constraints 

J^=nu, h + U^ = 0, (2.2) 

where u(x) is an m-component potential and h(x) is some m-component source term 
and 

nw^f^^y U^^-V-G + g (2.3) 



u 

serve to define the operators Fl and U. The differential constraints imply the key property 
that on any domain F 

f gj^+hu = [ {Gn)u, (2.4) 

Jr JdT 

as follows from integration by parts, where 

^■^^G-F + g-f, G-F = Tr(G^F), (2.5) 

and n is the outward unit normal to the surface dV. (Note that (12.51) is an inner product 
only when Q and JF are real, since it does not involve complex conjugation). Here, 
following standard notation, we define 

{Vu}„^|i, {V.G},^i:^. (2.6) 



Now suppose the fields T and Q are linked by the constitutive relation 

^(x) = Z(x)^(x), (2.7) 

where Z(x) = Z'(x) + iZ"(x) is a complex symmetric linear map (which is represented 
by a complex symmetric mid + 1) x m[d + 1) matrix when Q and T are represented by 
vnid + 1) component vectors). The map Z can be expressed in block form as 

K- m)- P'«) 

where (for each x) L(x) is a symmetric linear map on the space oi m ^ d dimensional 
matrices, K(x) is a linear map from m component vectors to m x c? dimensional matrices, 
K^(x) is a linear map from m x d dimensional matrices to m component vectors, and 
M(x) is a symmetric linear map on the space of m component vectors. Thus, written in 
terms of components, (12.71) implies 

Gij = LijkfFkf + Kijkfk, 
Qi = KkiiFki + Mikfk, (2.9) 

in which sums of repeated indices are assumed. Since L and M are symmetric maps, we 
have 

Lijke = Lkuj, Mik = Mki- (2.10) 

The constitutive law and differential constraints imply that u satisfies 

h + U(Znu) = 0, (2.11) 

or equivalently, 

V ■ (LVu + Ku) = h + K^Vu + Mu. (2.12) 

Our analysis thus applies to equations of this form, with L and M having the symmetries 

Let (u, JF, ^) be a solution of (12.21) and (12. 7p . Taking real and imaginary parts, the 
differential constraints imply 

j^' = nu', h' + u^' = o, j^" = nu", h" + u^" = o, (2.13) 

where the prime denotes the real part and the double prime the imaginary part, and the 
constitutive law can be written in the form 

?)-(i:-^z")(;')- '-' 

To make use of this representation let us further suppose that the imaginary part of 
Z satisfies 

Z"(x) > al, Vx e r, (2.15) 

7 



for some a > 0. Extending the ideas of ?) consider the real valued quadratic functional 

nu'\ /Z" TJ \ /nu' 



'2<-' -' = irlnu"j ■ (z' -Z"j inu"j +2(h".u +h'.u") 

= A(nu)-Znu + 2h-u]", (2.16) 

subject to the boundary conditions that u' = u' and u" = u" on dV . Let s(x) be a real 
valued m component vector field with s = on dV . Now Q(u'+s, u") is a convex quadratic 
function of s (with a linear component) while Q(u', u" + s) is a concave quadratic function 
of s (with a linear component). Therefore on the space of fields satisfying the boundary 
conditions the pair (u,', u") which are at the saddle point of the variational problem 



infsupQ(u',u"), (2.17) 



u' „// 



exists and is unique and corresponds to the stationary point of Q(u', u"), and is also at 
the saddle point of the variational problem 

supinfQ(u',u"). (2.18) 

u" u' 

Now let us prove we have the variational principle 

g(u',u") = inf supQ(u',u") = supinf Q(u',u"), (2.19) 



where the infimum and supremum are over fields with u' and u" with any fixed values 
Uq and Ug at the boundary dV and u' and u" are the solutions with these prescribed 
boundary conditions. 

Let (u, JF, ^) be a solution of (12.21) and (12. 7p . As s = on c^F, in a similar way as 
(12.41) was derived, we get 

I Q" ■ (ns) + h" ■ s = 0, I Q' ■ (ns) + h' ■ s = 0. (2.20) 

Thus we see that 

Q(u' + s, u") = g(u', u") + /" 2Q" ■ (ns) + 2h" ■ s + (Hs) ■ Z" n s 

= g(u', u") + /"(ns) • Z" n s, (2.21) 

and in the same way we get 

Q(u', u" + s) = Q(u', u") - /"(ns) • Z" n s. (2.22) 



So we see that the saddle point of f l2.17p is at (u', u") = (u', u"). 

Conversely, let (u', u") be a saddle point of (I2.17P satisfying some boundary condition 
u' = Uq and u" = Uq on dT. Consider 

infsupg(u' + s,u" + t), (2.23) 

s t 

where the infimum and supremum are taken over fields s and t satisfying s = t = on 
dr. A necessary condition for the saddle point to be at s = t = is that the first-order 
derivative of the functional Q vanishes: that is 

,' (nl) ■ iff) ^'''"■'^ '"'■'-'>' (2-24' 

for any s and t satisfying s = t = on SF, where we have introduced 



G"\ ^ f Z" Z' \ / nu' 

g' - [z' -z" J ■ V nu" 



(2.25) 



^'=(?.'V S"=(^'l]^ (2-26) 



Upon setting 



the condition (12.241) reduces to 

f G" ■ Vs + (g" + h") ■ s + G' • Vt + (g' + h') ■ t = 0. (2.27) 

Integrating this by parts gives 

/"(-V ■ G" + g" + h") • s + (-V ■ G' + g' + h') ■ t = 0, (2.28) 

and this will be satisfied for any s and t satisfying s = t = on SF if and only if 

h' + ug' = 0, h" + ug" = 0. (2.29) 

Thus the Euler-Lagrange equations associated with the saddle point variational principle 
coincide with the original equations (12. 2p . 

There is also a dual saddle point variational principle, analogous to the Thompson 
variational principle. By taking real and imaginary parts of the constitutive relation 
J-' = K.g where K = Z~^ has negative definite imaginary part we obtain 

-;)^(-_-)(aO- 

The dual variational principle involves the quadratic form 

fi(G'.G") = /(|:).(^: _^„)(|,)=/(g.K£)". (2.31) 

9 



where £' and Q^' are given by 

£'=(v.^-h')' ^"=(v-J-h")' ^^-^^^ 

and G' and G" satisfy the boundary conditions 

G'-n = G'-n and G" ■ n = G" ■ n on ^F, (2.33) 

where G' and G" are the real and imaginary parts of G. Notice that the form of Q_ given 
by fl2.32p ensures that h + U^ = 0. inside F. By similar argument to the one establishing 
(12.191) . we have the dual variational principle 

R{G\G") = supinf i?(G',G") = inf sup/?(G', G"), (2.34) 

G' ^" G" Q' 

where the supremum and infimum are over fields satisfying (12.331) . 

Again building upon the ideas of ?) let us rewrite the constitutive relation in the 
form 



where straightforward algebra shows that 



Z" + Z'(Z")-iZ' Z'(Z")-^ 

(z")-iz' (z")-^ 



(2.36) 



For all X G F the matrix C{'x) is positive definite since, using (12.141) and (12.351) . we see 
that the associated quadratic form 



= T' ■ TI'T' + T" ■ TI'T" (2.37) 

is non-negative, and zero only when T' = T" = 0, i.e. when J-"' = Q' = 0. 
Now consider the convex functional 

Y{u:, ^') = / ( % )-^{%)+ 2h" • u', (2.38) 



where Q' is given by 



^'=( V-^'-h'l' ^2.39) 



(to ensure that h' + U^' = inside F) and u' and G' are real and satisfy the boundary 
conditions that 

u' = u[), and G' • n = G[) • n on ^F. (2.40) 

10 



Let us prove we have the variational principle 

F(u',G') =infinfF(u',G'), (2.41) 

where the infimums are over all trial fields u' and G' satisfying the boundary conditions 
(12301) • Let (u, g) be a solution of dOD and (122D with u' = u[) and G' ■ n = Gg ■ n on 
dr. Suppose we are given any m component vector field s(x) with s = on SF and also 
a. m X d matrix valued field T(x) with T ■ n = on SF. Then (I2.20p and integration by 
parts implies that 



f 2g" ■ (ns) + 2h" ■ s - 2T ■ (nu") = O, where T= ( ^ J 



From this identity it follows that 



nu' + s,G' + T) = y(u',GO + ^(^^)-/:(^^) 



(2.42) 



(2.43) 



which implies (u', G') is a mimimizer of the variational principle. 

Conversely, let (u', G') be a mimimizer of the variational principle, satisfying some 
boundary conditions u' = Uq and G' ■ n = Gg • n on 9F. Consider 

infinfr(u' + s,G' + T), (2.44) 

s T 

where the infimum is over all fields s and T satisfying s = and T ■ n = on 5F. A 
necessary condition for the infimum to be at s = T = is that the first-order derivative 
of the functional Y vanish: that is 

,2(rT)-(j")+2h"-s = 0, (2.45) 

for any s and T satisfying s = and T ■ n = on dV, where we have introduced 



Upon setting 



r'=(]','X ^"=(5'')' (2-47) 



the condition (12.451) reduces to 

f G" ■ Vs + (g" + h") ■ s - T ■ F" - (V • T) ■ f" = 0. (2.48) 

Integrating this by parts gives 

/" (- V ■ G" + g" + h") ■ s - T ■ (F" - Vf ") = 0, (2.49) 

11 



and this will be satisfied for any s and T satisfying s = and T • n = on SF if and 
only if 

h" + ug" = 0, j^" = nf". (2.50) 

Upon identifying f" with u" we see that the Euler-Lagrange equations associated with 
the minimization variational principle coincide with the original equations (12.21) . 

The convex variational principle we have just stated is, as usual, associated with a 
dual variational principle: 

Y(u", G") = inf inf F(u", G"), (2.51) 



where 



F(u", G") = / ( £„ ) . £-1 ( £„ ) - 2h' . u", (2.52) 



in which 

£"=(v.#_h")' •'■''* 

and the infimum in (12.511) is over real fields u" and G" with u" = u" and G" ■ n = G" ■ n 
on dV. Since 

\-z'{z")-^ z" + z'(z")"^z7 ' ^ ' 

the expression for ^"(11", G") can be written equivalently as 

?(u", G") = ^ ( ]^" W ( ]^'' ) - 2h' • u". (2.55) 

These two minimization variational principles are not the only ones one can derive. 
Indeed, there is a continuous two-parameter family of related variational principles. The 
differential constraints (12.21) and the constitutive relation (12.71) imply 

^ = n5, h + U^ = 0, ^(x) = Z(x)^(x), (2.56) 

where 

S = e'"u, ^ = e*"^, G = e'("+^)^, h = e'^^^+^^h, Z(x) = e''^Z(x), (2.57) 

in which r and 6 are real parameters. So the fields with tildes on them satisfy the same 
equations as those without tildes, and the variational principles directly apply provided 
9 is chosen so that ^ 

Z"(x) > al, Vx G r, (2.58) 

for some a > 0. With the subsequent replacement (12.571) we obtain a new set of varia- 
tional principles parameterized by r and 6. Observe that even if (12.151) is not satisfied for 
any a > we still may be able to find a range of values of 6 for which (I2.58P holds, and 

12 



thus for which we can obtain variational principles. With the particular choice 6 = and 
e*"^ = —i the variational principle (I2.4ip gets mapped to the dual variational principle 

All the previous minimization principles need the fundamental coercivity assumption 
(12.151) . If Z"(x) = in some part \l/ of the domain F they cannot be applied directly. 
However we can still obtain variational principles by taking appropriate limits. If, in 
some region \I^ C F, Z"(x) is very small in the sense that 

61 > Z"(x), Vx G ^, (2.59) 

for some very small 6 > 0, then within this region 



^' ) ■ ^ ( _3' ) ~ [z' n u' - £'] • (Z")-1Z' n u' - £']. (2.60) 

If \1/ = F then the variational principle reduces to a variational principle similar to (I1.12p 
which is not particularly interesting. If \1/ 7^ F the variational principle will only be useful 
if one takes ^' ?« Z' Fl u' in \l/. The variational principles still hold in the limit as 5 — > 
provided one first takes 

£' = Z'nu', VxG^, (2.61) 

corresponding to taking an exact solution of (12. lip in this region (but with (u', G') not 
necessarily coinciding with the exact solution (u', G') associated with the given boundary 
conditions). Then as 5 ^ the expression for Y reduces to 

Outside \l/ we still require that Z"(x) > al for some a > 0. 

We now show how these variational principles may give useful information about the 
Dirichlet to Neumann map on the boundary dT. Let us assume h = and introduce the 
m-component complex vector fields 

qo = G(x) ■ n, Uo(x) = u(x) for x G dT. (2.63) 

The relation between qo and uq must be linear and we can write qo = Nuq which defines 
the Dirichlet to Neumann map N. If Uq and Uq denote two boundary conditions for 
u and Q^^\ J-'^^\ and qj/ , j = 1,2, denote the associated fields then it follows from the 
key property and the symmetry of Z(x) that 

JdV Jv JV 

= lQ'''-r^''=l q?)-uW = (u«,Nu?)), (2.64) 
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which imphes the map N is symmetric. It has positive semidefinite imaginary part since 

(u[„ N"u[,) + (u;;, N"u[;) = / q^ u[, - q[, ■ < 

Jar 

= I g" -T' -g' ■ T" = I T' ■ 71' T' + T" ■ TI'T" > 0. (2.65) 

Assuming N" is invertible, we can write the relation qo = Nuq as 



where 



^ _ /N" + N'(N")-iN' N'(N")-^ \ ., .7. 

^ -y (N")-iN' (N")-i yl ■ ^ ' 

Thus the map Af giving the boundary fields qg and Uq in terms of qg and Uq must have 
this special form. 

Now, using the key property, Y{u', G') can be expressed (when h = 0) in terms of 
the boundary fields: 

Y{u',G')= I g".:F'-g'.:F"= f q^ u;, - q;, ■ < 

Jr Jar 



Un \ . /• / Un 



and thus bounds on Y{u', G'), obtained from the variational principle, give upper bounds 
on the quadratic form associated with Af. Bounds on Y{u", G"), obtained from the dual 
variational principle, give lower bounds on the quadratic form associated with A/*. 

3 Application to Acoustics and Elastodynamics 

Let us begin by studying the acoustic equation in d dimensions: 

-V-p-^VP={l/K)uj^P, (3.1) 

where -P(x) is the complex pressure, /t(x) is the bulk modulus, p(x) is the density tensor, 
and uj is the (fixed) frequency of oscillation. [As is well known this equation with d = 2 is 
also applicable in cylindrical bodies, for antiplane elastodynamics where the displacement 
field is directed parallel to the cylinder axis, and for electromagnetism when the fields 
are transverse electric (TE), where the electric field is directed parallel to the cylinder 
axis, or transverse magnetic (TM), where the magnetic field is directed parallel to the 
cylinder axis]. Normally one expects k, to have a negative imaginary part (due to bulk 
viscosity), and p to be pi where p(x) is the positive mass density. However, viewed 

14 



as the effective density tensor of a composite material at a given frequency u;, p can be 
anisotropic (?) and, at least in the context of elastodynamics, can be complex valued with 
a positive semi definite imaginary part (?) with a real part that is not even necessarily 
positive [as established by ?), ?), and ?) in the context of antiplane elastodynamics, 
where 1/k plays the role of density, and by ?) and ?) in the context of three-dimensional 
elastodynamics]. Such unusual effective densities are due to the fact that different parts 
of the microstructure can experience different accelerations, and can be out of phase with 
the overall applied force if this relative motion dissipates energy. 

Given the complex pressure field -P(x) the associated complex velocity field of the 
fluid is 

V = -i{ujpY^WP. (3.2) 

[The associated physical pressure and physical velocity are respectively (e~*^*P(x))' and 
(e~*'^*v(x))' where the prime denotes the real part and t denotes the time.] Upon multi- 
plying ( 13. ip by e*^ and comparing it with ( 12. 12^ we can make the identifications 

u = P, L = -e^V"\ K = 0, M = e'^uj'^/n, h = 0, G = -ie'^uw. (3.3) 

Thus the potential u is a scalar, corresponding to the case m = 1. We choose 9 such that 
the imaginary part of 

e'^uj'^/K 

satisfies ( I2.15p . When p is real and positive definite, as is typically the case, then we 
need to choose 9 in the range Q > 9 > —it. If in addition k, is almost real and positive, 
with a small negative imaginary part, then we should choose 9 appropriately small and 
negative. 

Explicit expressions for the various variational principles in the acoustic case can of 
course be obtained by making the substitutions (13. 3p in the relevant equations in the 
preceding section. For example, if p has a positive definite imaginary part (which is 
more applicable in the mathematically analogous TE or TM electromagnetic problems) 
and fi:(x) has a negative imaginary part, then with 6* = we have the variational principle 

r(P',v") =infinfr(P',v"), (3.5) 

where 

v(P' ^r"\ - f I ' . .^, , , , . ,^ , 

-ujv" / I -cuv" \ ~V -w" I -V ■ V 



.e,,2/.. ) (3.4) 



n£', x") ^ / ( ^€. V ^ r ^€. V r ^^:.i ■ ^ i ^/.. i , (3.6) 



in which 

n 



r" + r'(r")-^r' r'(r")-^ 



f^_ik+[k)/kk/k. , . 
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and r = —p~^, and k = I/k. In this variational principle the boundary values of P' and 
v"-n are fixed, and the infimums are attained at fields P' and v" associated with solutions 
of the acoustic equation satisfying the prescribed boundary conditions. That the Euler- 
Lagrange equations of (13. 5p coincide with (13.11) and (13. 2p follows from the analysis of 
(I2.44I) - (I2.50I) . and can also be shown directly. 

When r" is very small within all of F then we have 



( ^S ) ■ ^ ( Ir ) " 1"''^'^ - "^"1 ■ ^'"r>'^^ - -^1. 



(3.i 



and the variational principles will only be useful if we take trial fields with uV' ~ r'VP'. 
The variational principles still hold in the limit as r" -^ provided we first choose 
v" = r'VP'/u;. Then, when p is real, we have the variational principle 

r(P') =infr(P') (3.9) 

where 

In this variational principle the boundary values of P' and n ■ p^^VP^ are fixed, where 
the latter corresponds to fixing the boundary value of v" ■ n. The infimum is attained 
at a field P' associated with solutions of the acoustic equation satisfying the prescribed 
boundary conditions. To show this, let P' be a minimizer and set P^ = P' + S where 
S = and n ■ p^^VS = on dT. A necessary condition for the minimum to be at S* = 
is that the first derivative of the functional Y vanish: that is 

,,J")-(v.p-'W)^°' <^-"' 

where we have introduced 

Integrating by parts the constraint (13.111) gives 

f S{ujU + V ■ p'^WP!') =0, (3. 13) 

and this will be satisfied for any S with S* = and n ■ p^^VS = on c?r if and only if 

f/ = -V ■ p- VP'Vcj. (3.14) 

A straightforward calculation (similar to the equivalence of (12. 7p and (I2.35P ) shows that 
(I3.12P is equivalent to (13. ip with this value of U. Thus the acoustic equation is the 
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Euler-Lagrange equation of (13.91) . when P" is defined by (13.121) . Tliis variational principle 
remains valid if we replace P by e^'^P where r can be any constant. 

In a similar way, if p has a positive definite imaginary part but k is real (which 
is more applicable to the analogous TM electromagnetic problem, where the magnetic 
permeability plays the role of 1/k) then we have the variational principle 



r(v") = infF(v"), (3.15) 

v" 

where 

in which the boundary values of v" • n and /tV ■ v" are fixed, where the latter corresponds 
to fixing the boundary value of P'. 

One has to shghtly modify the analysis to apply it to the elastodynamic equations, 

- V ■ (C Vu) = b + cj Vu, (3. 17) 

in which u(x) is the complex displacement field, b(x) is the complex body-force density 
p(x) is the complex density tensor, and C(x) is the complex elasticity tensor. Given 
the complex displacement field u(x), the complex stress t(x) and complex momentum 
density p(x) are given by 

T = CVu, p = —iujpvL. (3.18) 

[The associated physical displacement, body-force density, stress, and momentum density 
are respectively (e'^'^^u)', (e'^'^^b)', (e'^'^V)' and (e'^'^'p)'.] Upon multiplying fl3T7j) by 
e*^ and comparing it with (I2.12p we can make the identifications 

L = -e^^C, K = 0, M = e^Vp, h = e'^b, (3.19) 

and (13.181) implies 

G = -e^V, g = ie'^uJY). (3.20) 

The problem is that C acting upon any antisymmetric matrix is zero, and so L and hence 
Z are singular, which is a problem for computing the inverses in the equation (I2.36P for 
C This problem is rectified in the standard way, by replacing Fl everywhere with Fl 
defined by 

Then the fields T and Q have the form (12. ip . with F and G being symmetric d x d 
matrices: thus JF and Q can be represented by d{d + 3)/2 component vectors. We choose 
6 such that the imaginary part of 



-e'^C 
e'^u'p 
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id, ,2^ I (3.22) 



satisfies (12.151) on the space of tliese fields, and we compute the inverse of Z" on the space 
of these fields. When p is real and positive definite then we need to choose 9 in the range 
TT > ^ > 0. If in addition C is almost real with a positive definite real part, and with 
a small negative definite imaginary part, then we should choose 9 appropriately small 
and positive. With Fl replaced by Fl the analysis of equations fl2.23p - fl2.29p and fl2.44p - 
fl2.50p still holds, and we recover the elastodynamic equations from the Euler-Lagrange 
equations associated with the variational principles. 

In the special case ^ = the functionals Q, R and Y reduce to 

<9(u', u") = / (-e • Ce + cu^u ■ pu + 2b ■ u)", 

i?(T', r") = A-T ■ C-V + (V ■ T + b) ■ (cuV)~'(V • T + b)]", 

y((u',r'))= [ -( t]-c^ '' 



r 



-T ^ X — / 



+ l(b' + v%o/J^((b' + v%o/J+*"-^' ^'■'''> 



where 

and 

C 



e = [Vu + (Vu)']/2, (3.24) 



c" + cicy^c' c'(c")-' 



are respectively negative definite and positive definite matrices. It follows directly from 
the variational principles of ?) and ?), applied to the time harmonic case considered 
here, that the stationary points of the functionals Q and R correspond to solutions of 
the elastodynamics equations, in agreement with the results of the previous section. 
When p" is very small we have 

■^^ ^ ^ — I ^^ ' ' ."iJ^J I U' I V7 _'l/ .2 //\-lr ,2/, / 



(b' + V . r2)l. ; ^ ( (b' + V" ^')/. ) « l-VV+b'+V.rf(^V)- l-V'M'+b'H-V.^! 

(3.26) 
and the variational principles only will be useful if we take trial fields with cu^p'u' + b' + 
V ■ r' approximately zero. The variational principle involving the functional Y still holds 
in the limit as p" -^ provided we first choose 

u' = -(p')-i(b' + V-T0/cu2. (3.27) 

Then, when p is real, we have the variational principle 

F(T') = infF(T'), (3.28) 

T' 



where 

F(r') = ^ - ( _-^, ) ■ C ( _^^, ) + 2b" ■ u', (3.29) 

in which u' and e' are given by fl3.27p and fl3.24p . In this variational principle the trial 
field r' can be any real symmetric matrix valued field with prescribed values of r' ■ n 
and p~^V ■ r' at the boundary dT, where prescribing p~^V ■ r' at dT corresponds to 
prescribing u' according to (13.271) . The infimum is attained at a field r' associated with 
solutions of the elastodynamic equation satisfying the prescribed boundary conditions. 

4 The saddle point variational principles for electro- 
magnetism 

A comparison of the continuum elastodynamic equation (I3.17P with Maxwell's equations, 
VxE = zcjB, VxH = -icjD+j, D = £E, B = ^H, (4.1) 

which can be rewritten (?) in the form 

- V ■ (CVE) = icuj + Lu'^sE, (4.2) 

where now 

in which E is the electric field, j is the free current e the electric permittivity tensor, /x the 
magnetic permeability tensor, and Cijm = 1 (-1) if {hj,^ is an even (odd) permutation 
of (1, 2, 3) and is zero otherwise, suggests that the preceeding analysis should also extend 
to three-dimensional electromagnetism. Upon multiplying (14. 2 p by e*^ and comparing it 
with (I2.12p we can make the identifications 

u = E, L = -e'^C, K = 0, M = e'^tu^e, h = itue'^y (4.4) 

The problem is now that C acting upon any symmetric matrix is zero, and so L and 
hence Z are singular, which is a problem for computing the inverses in the equation (I2.36P 
for jC. This problem is rectified by replacing Fl everywhere with Fl defined by 

nE.(FE-(VEn/2y (4.5) 

Then the fields JF and Q have the form (12. ip . with F and G being antisymmetric 3x3 
matrices: thus JF and Q can be represented by 6 component vectors. We choose 9 such 
that the imaginary part of 

-e'^C 
e'^iu'^e 
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2. I (4.6) 



satisfies (12.151) on the space of tliese fields, and we compute the inverse of Z" on the space 
of these fields. 

Let us assume, for simplicity, that 

£"(x) > ail, /x"(x) > asl, Vx e T, (4.7) 

for some ai > and 0:2 > 0. (By multiplying (14. 2 p by e*^ and redefining e, n and j if 
necessary.) Then we can take 6 = and (12.151) holds for some a > 0. From (14.11) and 
(14.51) we see that 



F \ . -, „ i^ 



J- = nE = ( J. 1 , with F = -\ B, -5i I , (4.8) 



-^3 B2 
-B 
-B2 Bi 



and the constitutive law (12. 9p implies 

G 



-Hs H2 
-H 
-H2 Hi 



^ = ( ^2j3 ) with G = -iuj\ Hs -Hi | . (4.9) 



Thus the key property (12. 4p reduces to 

/ cu2(B-H + E-D) +ic<;j-E= / -icj(H x n) ■ E = / -icj(E x H) ■ n, (4.10) 
Jr Jar Jar 

and holds for all fields E, D, B,and H satisfying the differential constraints 

VxE = ic<;B, V x H = -icuD + j. (4.11) 

Now the first relation (|2:20|) will hold if 

/ (G" ■ n) ■ s = f -icoiU" X n) • s = 0, (4.12) 

Jar Jar 

and for this it suffices that the tangential component of s vanishes at the surface dT. 
Similarly the second relation in (I2.20p will hold if this condition is satisfied. Therefore 
the infimum and supremum in (I2.17P can be extended to all u' and u" with the same tan- 
gential components as, respectively, E' and E". Similarly in the minimization variational 
principle (I2.4ip . or (I2.5ip . we only need require that u', or u", have the same tangential 
components as E', or E". 

With n replaced by n the analysis of equations fl2:23|) - fl2:29|) and (^M>-(^M> still 
holds, and we recover Maxwell's equations from the Euler-Lagrange equations associated 
with the variational principles. If we introduce the tangential components 

E,|(x) = E-n(n-E), H,|(x) = H - n(n ■ H) (4.13) 
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of the fields E(x) and H(x) at the boundary dT, then 

qo(x) = G-n = -zcj(H|| X n). (4.14) 

When j = we can write qo = NEy, thus defining the map N which governs the 
electrodynamic response of the body at the frequency u. The analysis of fl2.63p - fl2.68p 
shows that the map N is symmetric with positive semidefinite imaginary part, and can 
be bounded (when j = 0) using the minimization variational principles. 

Equivalently we can find the variational principles directly from Maxwell's equations 
(14. II) . The key property (basically Poynting's theorem) (14.101) is obtained by integrating 
over r the identity 

V • (E X H) = H ■ (V X E) - E ■ (V X H) = iuj{U ■ B + E • D) - j ■ E. (4.15) 

Motivated by the form of (14. 8 p and (14. 9 p let us redefine the fields 

which are subject to the differential constraints that 

J^ = nE, iujj + Ug = 0, (4.17) 

where the operators Q and U are defined by 

nE=(J^^Y Ug = VxG + g (4.18) 

so that the differential constraints (14.111) and fl4.17p are equivalent. The key property 
(14.101) now takes the form 



g-J^ + iuji-E= (ExG)-n, (4.19) 

r Jdr 

where ^•JF = F-G + f-g. The constitutive law Q = ZjF holds with Z redefined as 
and this can be reexpressed in the form fl2.14p . with 



-^B"\ /^B'\ /^^H"\ f-uW 

E' J ' ^ =^^ =[e" )^ ^ = [u'B') ' ^ = [u'B" 

(4.21) 
where these fields satisfy the differential constraints 

j^' = nE', -uj}" + ug' = o, j^" = nE", uj}' + ug" = o, (4.22) 
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implying, as a corollary of f l4.19p . that 

g' -Qs- uj" -8 = 0, g" -Qs + ui' -8 = 0, (4.23) 

for all real valued vector fields s such that the tangential component of s vanishes at the 
surface dT. Introducing the functional 

QiS, E") -/(§')•( % % )(§')+ 2-0' • E' - J" ■ E") 

= \ [(fiE) ■ ZfiE + 2iu;j ■ E]", (4.24) 

we see from (14.231) that 

Q(E' + s, E") = Q(E', E") + /"(fis) ■ Z"fis, 

Q(E', E" + s) = g(E', E") - /"(fis) ■ Z"fis. (4.25) 

Consequently the pair of fields (E', E") = (E', E") is at the saddle point of (12.171) and of 
(12.181) . where the infimum and supremum are over all vector fields E' and E" with the 
same tangential components at the boundary dV as, respectively, E' and E". Due to the 
diagonal structure (14.201) of the matrix Z(x) the formula (I4.24p for Q(E', E") reduces to 

Q(E', E") = /"[(V X E) ■ m(V x E) + c^^e • eE + 22cjj ■ E]", (4.26) 

where m = —n^, like e, has positive definite imaginary part. Hence we have the saddle 
point variational principles 

g(E',E") = inf supQ(E',E") = supinf Q(E', E"), (4.27) 

E' E" E" E' 

where 



Q(E', E") = g".r + g'- r + 2a;(j' ■ E' - j" • E") 

cj'(H' ■ B" + H" • B' + E' ■ D" + E" ■ D') + 2cj(j' • E' - j" ■ E'). 

(4.28) 

When j = we can use (I4.15P to express (5(E',E") just in terms of the tangential 
components of the fields at the boundary: 

Q(E', E") =uj \ (E" X H") ■ n - (E' X H') ■ n when j = 0. (4.29) 

22 



The result that the stationary point of the variational principle f l4.27p corresponds to 
a solution of Maxwell's equations follows from a variational principle derived by ?). The 
Willis variational principle is expressed in terms of the electromagnetic potentials V(x) 
and A(x), which are associated with E and B through the relations 

E = -VV + iujA, B = VxA. (4.30) 

and, in the time harmonic case, involves (to within a proportionality constant) the func- 
tional 

Q(Z,A) =uj^ f {[-{VxA)-in{VxA)+{VV-iujA)-e{VV-iujA)+2iVV-i/uj-2A-i]e-'''y, 
Jar 

(4.31) 

where, as usual, the prime denotes the real part. His variational principle states that 

for any choice of t the functional Q{V_,A) applied to trial fields V_ and A having any 

given fixed values at the boundary dT, has a stationary point at the potentials V and 

A solving Maxwell's equations, and having the prescribed boundary values. Setting 

E = —'VV_ + iuA one sees by integrating by parts that Q{V_, A) is just a functional of E 

alone, and when t is chosen with e"*"^* = —i one has the identity 

Q(Z,A)=g(E',E"). (4.32) 

Furthermore the values of V_ and A on dT determine the tangential (but not normal) 
components of E' and E" on dT. 

There is also a dual variational principle. By taking real and imaginary parts of the 
constitutive relation JF = K.Q where K = Z~^ has negative definite imaginary part we 
obtain 02.231) . Let us introduce the functional 

fi(HMi") = I (1) (^: _'^„) (1) = /(g. K£)", (4.33) 

where Q^ and Q^' are given by 

2' = -(j"-VxH")' £"='^(-j'+"fxH')' <^-^^' 

(to ensure that — cuj" + UQ^ = and tuj' + 0£" = inside F) and H" and H' have the 
same tangential components as, respectively, H" and H'. We then have the saddle point 
variational principles 

i?(H', H") = sup inf i?(H', H") = inf sup i?(H', H"), (4.35) 

H" M' H' H// 

where the supremum and infimum are over fields with the required tangential components 
at the boundary. Due to the diagonal structure (I4.20p of Z(x) the formula for -R(H',H") 
reduces to 

^(H', H") = / [(j - V X H) • e(j - V X H) + cj^H^xH]", (4.36) 
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where e = —s~^, like /j,, has positive definite imaginary part. From this expression one 
can see that -R(H', H") is a special case of the functional associated with the ?) dual 
variational principle. We also have that 

R{u', H") = I g" -r + g' ■ r = f uj\ii' ■ b" + h" ■ b' + e' ■ d" + e" ■ d'), (4.37) 

which when j = can be expressed in terms of the tangential components of the fields 
at the boundary: 

i?(H', H") =uj (E" X H") • n - (E' X H') • n when j = 0. (4.38) 

Jar 

Consider a medium which is locally isotropic and such that £(x) is purely imaginary 

(which may be a good approximation at low frequencies in conductive media) while /x(x) 

is purely real (and so (14.71) only holds in the limiting sense with ^2 = 0). We can write 

£ = il/{ujp), e = —iujpl, n = /il, (4.39) 

where the resistivity p(x) and the permeability /i(x) are purely real and positive. Then 
the expression ( 14.36^ reduces to 



^(H', H") = f ujpili' - V X Hf - Ij" - V X H'f ) + 2c^VH' ■ H" 



(4.40) 



and then when j = the saddle point variational principle (I4.35P corresponds to the ?) 
variational principle. (Note that j has a different meaning in that paper.) Now Maxwell's 
equations still hold if we make the replacements 

£^-ie, M^i/x, E^iE, B ^ iB, D ^ D, H ^ H, j^j, (4.41) 

(which is a special case of the transformation (I2.57p ) and with these replacements applied 
to (14.391) the expression (I4.36p reduces to 



^(H', H") = f -2ujp{}' - V X H') ■ (j" - V X H") + cu^dHf - |H'f ) 



(4.42) 



and when j = the saddle point variational principle (I4.35P then corresponds to the other 
Borcea variational principle. Multiplying e and dividing fj, by e*^, where > 6 > —tt/2 
leads to a continuous family of saddle point variational principles which interpolate the 
two Borcea variational principles. 

5 Minimization variational principles for electromag- 
netism 

To obtain a minimization principle we rewrite the constitutive law in the form (I2.35P 
where C is given by (12.360 . Then we introduce the quadratic functional 

r (E', H") = ^ (^ ^|! ) ■ £ (^ ^|! ) + 2ui' . E', (5.1) 
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where Q' is given by 



Q! = ^{ j"_ VxH" ^ ' ^^"^^ 



(to ensure that —coj" + UQ^ = inside F) and E' and H" have the same tangential 
components at the boundary dT as, respectively, E' and ujH". Given two vector fields 
s and t with zero tangential components at the boundary dT then (I4.23P and the key 
property fl4.19p imply 

f G" ■ (fis) + cjj' • s - T ■ (fiu") = 0, where T = uj ( _^ ^ ^ j . (5.3) 

As a result we deduce that 

r(E' + s,H" + t)=y(E',H") + ^(^^V£(^^V (5.4) 

Consequently, because C is positive definite under the assumption 04. 7p . we the mini- 
mization variational principle 

F(E',H") =infinfy(E',H"), (5.5) 

E' H" 

where the infimums are over vector fields E' and H" with the same tangential components 
at the boundary dV as, respectively, E' and H". Using (12.371) the minimum value of 
F(E',H") is 

F(E', H") = I uj\W ■ B" - H" ■ B' + E' • D" - E" • D') + 2cuj' • E'. (5.6) 

Now the power dissipated into heat in the body F, averaged over a cycle of oscillation is 

Vr(E,H) 



r L 



/ /Ee-^"*)' ■ ^^ '- + (He-'"*)' • ^^ '- + (Ee"*"*)' ■ (je"'"*)' dt 



27r7o dt ' ' dt 



- /cu(H' ■ B" - H" • B' + E' ■ D" - E" ■ D') + E' ■ j' + E" ■ j", 
2 Jr 



(5.7) 



which is proportional to (15.61) when j = 0. So the minimum value of Y has a physical 
interpretation when j = 0. 

Due to the diagonal structure (I4.20p of Z the formula (15. ip reduces to 

''®'s")-/(vxf'-r)-<vxf'-r 
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where ^(x) and A4(x) are the positive definite matrices 

■^- [ (m")-im' (m")-^ J ' ^ ' 

in which m = —f-i^. Also the constitutive law ( 12.351) reduces to 

So given any arbitrary prescribed values for the tangential components E' and H" at 
the boundary dT the minimum in the variational principle (15. 5p will satisfy the Euler- 
Lagrange equation and therefore will correspond to a solution of Maxwell equations. 
From the minimizing fields E' and H" one can, according to (I5.10p . determine 

E" = (s")-'s'E' + (£")"'(V X H" - j")/cu, 

H' = m'{m")-^¥L" - [m" + m'(m")"'m'](V x E')/cJ, (5.11) 

and in particular one can determine the tangential components of these fields at the 
boundary. 

It is often the case that n is real, while £ has a positive definite imaginary part. In 
this event we redefine Z, Q and j by multiplying them by e*^. Again there is a continuous 
two-parameter family of related variational principles. We rewrite Maxwell's equations 
(14.11) in the form 

VxE = zcjB, VxH = -zcjD+J, D = £E, B = /iH, (5.12) 



where 



E = e"*^*«E, D = e*(^~'^*°)D, H = e''(^~'^*°)H, B = e"''^*oB 



j = e^^^--*«^j, e = e''e, m = e'^, (5.13) 

in which the real constant to (which may be viewed as a time interval) can be chosen 
freely, while 6 must be chosen so that 

i"(x) > ail, /i"(x) > aal, Vx G T, (5.14) 

for some cti > and 02 > 0. Then the fields with the tildes on them also satisfy Maxwell's 
equations and the variational principles directly apply to them. With the subsequent 
replacement (I5.13P we obtain a new set of variational principles parameterized by 6 and 
to- When 6 is non-zero the minimum value of the minimization principles no longer has 
a physical interpretation as the time averaged power dissipation in the body. 
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In particular by choosing 6 = and to with e *'^*° = —i we recover the dual variational 
principle that 

is minimized at E" = E" and H' = H', under the constraint that E" and H' have the 
same tangential values as, respectively, E" and H'. A corollary is that 

W{E,H) = [F(E',H") + r(E",H')]/(4cu) (5.16) 

is minimized at E = E' + iE" = E and H = H' + iH" = H, under the constraint that 
E and H have the same tangential values as, respectively E and H. By making the 
substitutions fl5.13p in fl5.6p it is easy to check that the minimum value of W^(E, H) is the 
average power dissipation ly (E, H) given by (15.71) . even when j is non-zero. An appeal- 
ing feature of the expression (4.39) is that it remains invariant under the replacement 
E -^ e~*'^*°E, H — >■ e"""^*"!! for any real value of to. However a disadvantage is that 
the boundary conditions are overprescribed: if one imposes boundary conditions on the 
tangential values of E and H that are not compatible with the moduli £(x) and /ii(x) 
then the fields which minimize W^(E, H) will not satisfy Maxwell's equations. 
When ^", or equivalently m", is very small we have 

{"^-uS) ■ -^ ("^-uS) "^ f™'^ "" - " ""-^ ■ (™")"'KV X E' - cuH"] (5.17) 

and the variational principle is only useful if we take trial fields with H" ^ m'V x E'/u;. 
If we take H" = m'V x E'/ci; we still obtain a useful variational principle in the limit 
when (J, is real. From (15. 5p and (15.80 we have 

r(E') = infr(E'), (5.18) 

E' 

where 

^^-^ ^ X V-V X M-^(v"x EO/c^- rj '^ (-V X t^-\V~x E')/^- rj +2'^J'-^'- 

(5.19) 
In this variational principle the infimum is over vector fields E' with prescribed tangential 
components of E' and i-i~^V x E' at the boundary dT. The infimum is attained at a field 
E' associated with the solution of Maxwell's equations satisfying the prescribed boundary 
conditions. 
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6 Application to Electromagnetic Tomography 

Here we consider the tomography problem: how can one recover information about the 
functions £:(x) and A^(x) given measurements of the tangential components E||(x) and 
H||(x) of the fields E(x) and H(x) at the boundary SF? When it is known that the body 
r consists of small inclusions in a matrix this question has been answered by ?) : see also 
the review by ?) . For other geometries a different approach is needed and the variational 
principles provide this. First observe from fl4.15p and (15 .Tp that the power dissipation 
Vr(E, H), only depends on the tangential values of these fields at the boundary 



W{E, H) = -i / (Ejj X Hj|) ■ n + (Ej[ x Hj[) ■ n, 



(6.1) 



and, as could be expected physically, is the time-averaged value of the flux of the Foynting 
vector into the body F. Thus when j = 0, and trial fields E' and H" are chosen having 
the same tangential values at the boundary as, respectively E' and H", the inequality 

W{E, H) < r (E', H")/(2cu) (6.2) 

provides for a given £(x) and At(x) a bound on the possible tangential components of 
the fields E" and H'. Conversely if these tangential components have been measured the 
inequality provides one constraint on the possible values of £(x) and At(x). Additional 
inequalities on £(x) and ^t(x) may be obtained by chosing different trial fields, by working 
with other minimization variational principles in the family parameterized by tg and 6, 
and by conducting a series of physical experiments with various different sets of values 
of the tangential components of the fields E(x) and H(x) at the boundary dT. It is not 
necessary that the moduli £(x) and //(x) remain the same for each physical experiment 
provided we have some model for how they vary from experiment to experiment. For 
example, the frequency u could be different for the different experiments if we know 
that the permittivity has the form £(x) = Sq + icTo/uo (or some other known parametric 
dependence on uo). The objective is then to constrain the possible values of the real 
functions £o(x) and cro(x). 

One criticism of this approach is that it only utilizes information (in the case 6* = 0) 
about the time- averaged power dissipation in the body F. When t = 9 = Q the measured 
values of the functions E" and H' around the boundary dV are not used except to calculate 
the single scalar quantity iy(E, H). This seems like a tremendous waste of information. 
However some extra information from the measured fields can be incorporated if one has 
measurements of the boundary value fields Ey (x) and H|| (x) for n experiments indexed 
by the integer j = 1,2, ... ,n, with the moduli £(x) and ^J,{x.) remaining the same for 
each experiment. Then one can infer the boundary values Ey and Hy associated with a 
linear combination of these fields: 

n n 

E|| = 5^A,e[|^-\ H|| = 5^A,h[|^-), (6.3) 
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where Ai, A2, . . . A„ are a set of complex constants. Letting E and H denote the corre- 
sponding electric and magnetic fields inside F, it follows that H^(E, H) can be expressed 
as the sum 

n n 

iy(E, H) = 5^ 5^(a;.a', + \",K)w,, + (a;a', - K\'^)s,u 
j=i fc=i 

^ n n 

E E^^^-^'fc + KK)iW,, + w,,) + (a;a', - K\'^){s,, - s,,) (6.4) 
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/(Ej 


(i) 


xHj 


Kfc) 



n + (Ej[^^'^ X Hjf "^V n, 

n-(Ej[^^'^ xHJi^'^Vn. (6.5) 

which can be extracted from the experimental measurements of the boundary fields. Thus 
iy(E,H) depends not only on the n diagonal elements, Wjj, which physically represent 
the time-averaged power dissipation in each experiment, but also on the n(n — 1) elements 
Wjk + Wkj and Sjk — Skj, with j > k, which have no such physical interpretation. It 
makes sense to choose trial fields which are also linear combinations: 

n n 

E' = J2(x0'^y = E ^jE'^'^ - y-w^'\ 



n n 



H" = 5^(A,H(^))" = J2 A;H"(^) + A;'H'(^1 (6.6) 

i=i i=i 

Then y(E', H")/(2u;) — W^(E, H) must be a positive quadratic form in the 2n real vari- 
ables X'l, A'/, A2, A2 . . . A^, A", and the positive semi-definiteness of the associated 2n by 
2n matrix will give constaints on the functions £(x) and /x(x). 

We can also use the variational principle (I5.18P for tomography purposes. When /j, is 
real and j = the inequality 

Vr(E,H) <F(E')/(2cj) (6.7) 

applies. Often one knows in advance that the material inside F is non-magnetic, so that 
H = I inside F. Then from knowledge of PF (E, H) (16.71) provides a constraint on the 
possible values of the electrical permittivity £ for every suitable choice of the trial field E'. 
On the other hand, suppose the medium inside F was two phase with an unknown phase 
boundary and with two different values of fj, inside F. Then the variational principle (16. 7p 
will be useless unless the choice of trial field is carefully correlated with the unknown 
phase boundary, so that the tangential component of /x^^(V x E') is continuous across 
the phase boundary. 
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When j 7^ then we have the inequahty 

iy(E,H) < iy(E,H), (6.8) 

where Vr(E, H) is given by f l5.16p and the trial fields E and H are required to have the 
same tangential values as, respectively, E and H. Given measurements of the tangential 
values of E and H, and trial fields E and H we expect that this inequality should provide 
a useful constraint on the possible values of £(x) and /^(x). 
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